Polaronic and nonadiabatic phase diagram from anomalous isotope effects 
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Isotope effects (IEs) are powerful tool to probe directly the dependence of many physical properties 
on the lattice dynamics. In this paper we invenstigate the onset of anomalous IEs in the spinless 
Holstein model by employing the dynamical mean field theory. We show that the isotope coefficients 
of the electron effective mass and of the dressed phonon frequency are sizeable also far away from the 
strong coupling polaronic crossover and mark the importance of nonadiabatic lattice fluctuations in 
the weak to moderate coupling region. We characterize the polaronic regime by the appearence of 
huge IEs. We draw a nonadiabatic phase diagram in which we identify a novel crossover, not related 
to polaronic features, where the IEs attain their largest anomalies. 

PACS numbers: 71.38.-k,63.20.Kr,71.38.Cn,71.38.Ht 



The relevant role of electron-phonon (el-ph) interac- 
tions on the properties of complex materials as the high- 
T c superconductors has been recently revived by a num- 
ber of experiments. In particular, the finite, and yet un- 
explained, anomalous IEs on the penetration depth0], 
on the pseudogap temperature[2j , and on the angle- 
resolved photoemission spectra|2j in high-T c supercon- 
ductors open a new challenge in the understanding the 
role of the el-ph interaction in these materials. Sizable 
isotope effects on magnetic and charge-ordering critical 
temperatures in manganites point out the relevance of 
the el-ph coupling also in these materials 4] . 

Despite that the el-ph problem has been thoroughly 
studied and that many of its manifestations are now well 
understood, only few partial studies have been devoted to 
IEs and to their significance in relation to the underlying 
nature of the el-ph interaction [HQ Yet, the prediction 
and observation of IEs on different physical properties 
represent a powerful tool to assess the role of the el-ph 
interaction in many materials. For instance, the finite 
isotope shift on the penetration depth [jj, and hence in- 
directly on the effective electron mass m*, is of partic- 
ular interest since it contrasts the conventional Migdal- 
Eliashberg (ME) scenario which predicts strictly zero IE 
on this quantity. In this perspective the understanding 
of finite IEs on m* cannot rely on the ME framework and 
more general approaches are then required. 

A modern tool of investigation which overcomes the 
limitations of ME theory is the so-called dynamical 
mean field theory (DMFT), a nonperturbative method 
which neglects spatial correlations in order to fully ac- 
count for local quantum dynamics, and becomes exact 
in the infinite coordination limit [3- In the case of el- 
ph interactions, this approach allows us to study with 
equal accuracy all coupling regimes, and to fully in- 
clude phonon quantum fluctuations which are only par- 
tially taken into account in ME approach according to 
Migdal's theorem. DMFT has been successfully em- 



ployed in the study of multiphonon effects, polaron in- 
stabilities, metal-insulator transitions (MIT) as well as 
quasiparticle regimes in the Holstein electron-phonon 

systemHUEtlll 

Considering the Holstein model as a paradigmatic lat- 
tice model for el-ph interaction, we discuss the anomalies 
of the IEs on electronic and phononic properties arising 
in purely el-ph systems. We employ DMFT to span the 
whole phase diagram determined by the dimcnsionlcss 
el-ph coupling A and by 7 = Wo/t, the "adiabatic" ratio 
between the phonon energy luq and the electron hopping 
rate t. Increasing A gradually induces a crossover from 
free carriers to polarons with reduced mobility, while 7 
measures the relevance of quantum lattice fluctuations. 
Particular attention will be paid on the dependence of 
the IEs on the quantum lattice fluctuations triggered by 
the finite adiabatic parameter 7 7^ 0. 

We show that a sizable negative isotope coefficient on 
the electronic mass m* characterizes the whole param- 
eter range, while a divergence of a m * is recovered only 
approaching the MIT polaron transition at A c and 7 = 0. 
Similar features are found for the renormalized phonon 
frequency Qq which displays an isotope coefficient ao 
significantly different from the ME limit ct<n = 1/2. 
Based on the dependence of the isotope coefficients on the 
adiabatic ratio 7, we draw a phase diagram wherein we 
identify, beside the strong-coupling polaronic regime 10],: 
i) a nonadiabatic perturbative regime where IEs increase 
with 7, and ii) a complex nonadiabatic regime where the 
anomalies of the IEs decrease with 7 and approach the 
Lang-Firsov predictions in the too/t — ► 00 limit. 

In this paper we are interested in the continuous evo- 
lution of the el-ph properties from the quasiparticle to 
the polaronic regime. It is well known that the Holstein 
model undergoes various instabilities leading to super- 
conductivity, charge-density-wave ordering and bipolaron 
for mat ion [12]. In order to focus on the metallic proper- 
ties and to clarify the origin of anomalous IEs in this 
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regime, we consider here a half-filled spinless Holstein 
model, which enforces the metallic character in the whole 
A-7 space (for 7 7^ 0) [T^j. Our Hamiltonian reads: 
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where cj (cj), aj (ai) are creation (annihilation) operators 
for electrons and phonons on site i, respectively, rij = c\ci 
is the electron density and g is the el-ph matrix element. 

The physical properties of Q ar e governed by two mi- 
croscopic parameters: the el-ph coupling A = 2g 2 /u>ot 
and the adiabatic ratio 7 = wo/t. The standard Landau 
Fermi-liquid (FL) picture is sustained by the ME theory 
in the adiabatic limit 7 = for A < A c , while for A > A c 
the FL regime is destroyed due to the polaron localiza- 
tion. For 7 > 0, this sharp transition becomes a smooth 
crossover between well and poorly defined quasiparticle 
properties. This rich phenomenology is reflected in the 
appearances of anomalous IEs on various quantities. Let 
us consider for example the isotope coefficient on the 
effective electron mass m*: a m - = — dln(m*)/dln(M), 
where M is the ionic mass. Since g oc 1/yMwo and 
loq oc 1/Vm, the el-ph coupling A is independent of M 
and a m * can be rewritten as: 



1 dhi(m* / m) 

2 dln( 7 ) 



(2) 



According to the FL description, m* can be expressed in 
terms of a mass-enhancement factor f m * : 



m*/m = 1 + / m »(A,7). 



(3) 



In the adiabatic regime 7 -C 1, fm* can be expanded in 
powers of 7, m*/m~ 1 + f m * (A, 0) + 7/^1* (A), and hence 



7 m 
2 m* 



(4) 



The isotope coefficient thus increases with 7 and it cor- 
rectly reproduces the ME result a m * = when 7 — > 0. 
Such an increase is indeed found by calculations based 
on a perturbative expansion in 7 [a, 0], whose validity 
is however limited only to weak A values. In the oppo- 
site antiadiabatic limit 7 — > 00 the Holstcin-Lang-Firsov 
approximation gives m* /m ~ exp(A/27), leading to 



A_ 

'47' 



(5) 



which decreases as 7 gets higher, as opposed to to the pre- 
vious case of Eq.(£j}. The two limiting cases for 7 <C 1 and 
7>1 suggest that for fixed A < A c the strongest isotope 
shifts lie in the intermediate nonadiabatic region 7 < 1, 
which can be investigated only by non-perturbative tools 
as the DMFT approach. 

In this work we consider JJ| on an infinite coordina- 
tion Bethe lattice and use exact diagonalization (ED) to 
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FIG. 1: Effective electron mass m*/m (a) and renormalized 
phonon frequency Qo /^>o (c) as function of A for 7 = 0.1 (solid 
line), 7=1 (dashed line) and 7 = 10 (dot-dashed line). The 
corresponding isotope coefficients are shown in panels b,d. 



solve the impuri ty p roblem that DMFT associates to the 
lattice modelQ, Il3j ■ As customary, the Anderson model 
is truncated by considering N s impurity levels, and a 
cutoff on the phonon number is imposed on the infinite 
phonon Hilbert space. The DMFT self-consistency is im- 
plemented in the Matsubara frequencies u> n = (2n+l)nT 
where T is a fictitious temperature. The evaluation of IEs 
is a particularly difficult task since it requires a high ac- 
curacy on both the electron and phonon properties and of 
their dependence on 7. In particular, a number of phonon 
states up to 100 and T/t as small as 1/1600 were needed 
to ensure reliable and robust results. The number of im- 
purity levels has been fixed at N s — 9, having checked 
that no significant change occurred for larger N s . We 
compute the electron self-energy E(a>„) and the phonon 
Green's function D(oj m ) which yield the effective electron 
mass m* I'm = 1 — S(w n= o) / 'ttT and the renormalized 
phonon frequency Qq as (Qq/uq) 2 = —2D~ 1 (u> m= o)/t*Jo. 
The corresponding isotope coefficients are obtained by 
means of a finite shift A7/7 = 0.15. 

In Fig. ^ we show m*/m and flo/^o and their corre- 
sponding isotope coefficients as a function of the el-ph 
coupling A for 7 = 0.1, 1, 10 which are representative re- 
spectively of the quasi-adiabatic, nonadiabatic and an- 
tiadiabatic regimes. The polaron crossover is reflected 
in a strong enhancement of m*/m as A increases. The 
crossover occurs at A = 1.18 for 7 = 0, and moves to 
larger coupling as 7 increases, until the antiadiabatic 
regime, in which the crossover roughly occurs when the 
average number of phonons a 2 = A/27 > 10, Hoi IT^ |. 
As a consequence, at fixed A the effective mass becomes 



FIG. 2: (color online) Phonon spectral function (solid lines) 
for different values of the el-ph coupling A from to 10 with 
steps of 0.2. The thick red line represent the renormalized 
phonon frequency flo and the thin dashed line u>o. 

smaller as 7 increases, implying a negative isotope coeffi- 
cient a m » as reported in Fig^b). The polaron regime is 
thus identified by the huge negative values of a m * . Just 
as in the case of m*/m, increasing 7 makes the variation 
of a m * smoother and shifts it to larger A. 

A similar behavior is found for the renormalized 
phonon frequency flo, as shown in Fig^c), where the 
polaron instability is reflected in a sharp phonon soften- 
ing in the quasi-adiabatic case (7 = 0.1) close to polaron 
crossover. Once again, the softening is weaker as 7 gets 
larger. For fixed A this leads to an anomalous isotope 
coefficient an Q > 1/2, as reported in Fig^d). 

Some words are worth to be spent about the appar- 
ent phonon hardening accompanied by the corresponding 
decreasing of an a as reported in FigsQIc,d) f° r large A. 
This anomalous phonon feature is just a consequence of 
our "static" definition of renormalized phonon frequency 
{yio/ujo) 2 — —2D~ 1 (oj m =o)/uJo. This definition corre- 
sponds to describe the full phonon spectrum as a single 
(^-function at frequency f2g = — l/(iru)o) J du>Im.D(uj + 
i0 + )/uj. Such description becomes less representative as 
the spectral function acquires a complex structure. 

In Fig|3we report the evolution of the phonon spectral 
function by increasing the el-ph coupling A for 7=1. The 
real frequency phonon propagator is directly computed 
in the ED scheme. Note that only gross features can 
be extracted because of the discreteness of the impurity 
model. The average phonon softening (thick red line) at 
small A stems from a transfer of spectral weight from lu ~ 
ujq to a low energy peak which exponentially approaches 
w = 0. An opposite behavior occurs for strong el-ph 
coupling where the lattice potential is a double-well with 
energy barrier AE ^> luq. This gives rise to a second 
peak at frequency luq yielding a hardening of the averaged 
phonon frequency Similar considerations lead to the 
decrease at large A of the IE an . 

The study of the A dependence of m*/m, tto and their 
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FIG. 3: Dependence of a m * and an on the nonadiabatic 
parameter 7. Curves are plotted for the el-ph coupling A = 
0.6,0.8, 1.0, 1.2, 1.4, 1.6 from top to bottom (bottom to top) 
in left (right) panel. Filled black circles mark the adiabatic 
ME limit. 



respective IEs highlights in the clearest way the appear- 
ance of giant IEs close to the polaronic crossover. How- 
ever, interesting anomalies and sizeable deviations of IEs 
from ME predictions appear far from polaronic regime in 
a region where the system preserves good metallic prop- 
erties. These anomalous features can thus be connected 
to the onset of nonadiabatic effects in the el-ph properties 
in the spirit of a g eneralization of the FL liquid picture 
as in Eq.J3J) 5. Il5l ITtij. To this aim we analyze in more 
details the dependence of the IEs on the quantum lat- 
tice fluctuations triggered by the finite adiabatic param- 
eter 7^0. In Fig. |21 we plot the 7-dependence of a m * 
and Qfj for different el-ph coupling ranging from weak 
(A = 0.6) to intermediate-strong coupling (A = 1.6). Two 
qualitatively different behaviors marked by the solid and 
dashed lines are identified. Dashed lines represent strong 
coupling A values for which the system undergoes a pola- 
ronic MIT for 7 — > signaled by the divergence of the iso- 
tope coefficients. Solid lines are representative of weak- 
to-moderate A values for which the system maintains its 
FL metallic character with isotope coefficients recover- 
ing their standard ME values a m * — and ao = 1/2 
as 7 — > 0. Note that for A = 1.2 we find that a m * — * 
and ev£i — ► 00 for 7 — * 0, in agreement with Ref. [llT ] 
which predicts two different critical values A^ h = 1.18 
and Xf = 1.328 where a double- well lattice potential and 
an electronic self-trapping respectively occur. 

The different physics underlying the regimes sketched 
by dashed and solid lines is reflected in the overall 7 de- 
pendence of the IEs. In the large coupling regime A > A c 
the isotope coefficients diverge for 7 — > and approach 
monotonically zero in the antiadiabatic limit. For A < A c 
on the other hand we find a more complex behavior with 
an initial increase of the anomalies of a m * and an as 
function of 7 up to a certain value 7*, above which the 
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FIG. 4: (color online) Nonadiabatic phase diagram defined by 
anomalous IEs. Filled red circles and black crosses represent 
respectively the minima and the maxima in the 7 behavior of 
cim* , Qn • Error bars are not related to the accuracy of the 
isotope coefficients but more to the flatness of a m * , &n as 
function of 7. The curve A c marks the polaron crossover as 
defined in the text. 

anomalies decrease and disappear in the antiadiabatic 
7 — > 00 limit in agreement with the Lang-Firsov predic- 
tions. In the nonadiabatic crossover region (A < A c and 
7^7*) significant anomalies of the IEs a m » ~ —0.3 and 
a n ~ 0.55 are observed even in the metallic regime far 
from the polaronic instability. It is worth to stress that 
the values of 7* strongly depends on the el-ph coupling 
A and that in principle a nonadiabatic crossover param- 
eter can be defined for both the electronic and phononic 
properties. 

As discussed above, the curve 7* separates a region in 
which a perturbative theory based on the adiabatic FL 
picture can be safely employed (7 < 7*) from a region 
where the Lang-Firsov approach is a more appropriate 
starting point for a I/7 expansion (7 > 7*). However 
around A c at 7 ~ both ME and Lang-Firsov approaches 
are expected to fails due to strong entanglement of elec- 
tron and phonon degrees of freedom || . 

The complex phenomenology of the IEs in the whole 
A-7 phase diagram is summarized in Fig. where the 
nonadiabatic crossover 7* is defined respectively by the 
minima (filled circles) and maxima (crosses) of the a m * , 
an as functions of 7. We note that, apart from a very 
small region close to the adiabatic limit, 7* lies on an 
universal curve for both the electron and phonon prop- 
erties reflecting the fact that electron and lattice degrees 
of freedom are strictly mixed in the nonadiabatic regime. 
This is indeed not true for 7 — > where 7* disappear at 
the two different values A£ h , Ajf above discussed. It is also 
worth to stress that the nonadiabatic crossover, pointed 
out by the 7 dependence of the anomalous isotope effects, 
is not related to the onset of polaronic effects, as shown 
in Fig. 01 where the polaron crossover at finite 7 is here 



pinpointed by the appearance of a bimodal structure in 
the lattice probability distribution function. For A > A c 
the polaron crossover is reflected in a strong degrada- 
tion of the metallic properties with an almost vanishing 
coherent spectral weight. The behaviors of A c and 7* 
define three regimes which can be roughly described as: 
I) a polaron region (brown dashed area in Fig. where 
electrons are almost trapped leading to significant lat- 
tice distortions. This effect is strongly removed by the 
lattice quantum fluctuation triggered by a finite loq lead- 
ing to giant IEs; II) a highly nonadiabatic region (white 
area) where the system is qualitatively described in terms 
of itinerant quasi-particle carrying along its hugely fluc- 
tuating phonon cloud. In this regime quantum lattice 
fluctuations relax the mixing between lattice and elec- 
tronic degrees of freedom leading to a reduction of the 
anomalous IEs; III) a weakly nonadiabatic region (grey 
filled area) where anomalous IEs are tuned by the open- 
ing of nonadiabatic channels in the el-ph interaction. In 
this region DMFT qualitatively confirms the results of 
the nonadiabatic theory described by vertex diagrams in 
a perturbative approach [5! Ho]]. 

In conclusion we have defined a phase diagram of the 
spinless Holstein model based on the anomalous phe- 
nomenology of the IEs on both electronic and phonon 
properties. In the metallic regime we identified a new 
nonadiabatic crossover 7* (not related to polaronic in- 
stability) between different weakly and highly nonadia- 
batic FL quasi-particle picture. The largest anomalies 
in the isotope coefficients are found in this intermediate 
crossover region. Experimental investigations of anoma- 
lous IEs on different physical properties represent thus a 
powerful tool to probe the complex nature of the el-ph 
interaction in real materials. Theoretical studies based 
on Density Functional Theory predict A 3 Cqq fullerides, 
cuprates and MgB2 to be roughly in the nonadiabatic 
crossover where the anomalies in the IEs are the most 
sizable. A rigorous generalization of our results in the 
presence of electronic correlation and away from the half- 
filling case is of course needed for quantitative analysis in 
these complex materials. Preliminary results show that 
a m * can be significantly suppressed in low filled systems 
even in highly nonadiabatic regime, in agreement with 
recent experimental results in MgB2 01 ■ 
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